Phase diagram of a polarized Fermi gas across a Feshbach resonance in a potential trap 
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We map out the detailed phase diagram of a trapped uhracold Fermi gas with population imbal- 
ance across a wide Feshbach resonance. We show that under the local density approximation, the 
properties of the atoms in any (anisotropic) harmonic traps are universally characterized by three 
dimensionless parameters: the normalized temperature, the dimensionless interaction strength, and 
the population imbalance. We then discuss the possible quantum phases in the trap, and quantita- 
tively characterize their phase boundaries in various typical parameter regions. 



I. INTRODUCTION 

Recent experiments with ultracold Fermi gases near 
a Feshbach resonance, through which the inter-atomic 
interaction can be tuned by varying the external mag- 
netic field, have attracted considerable attention [J . The 
latest advance in this direction features two very recent 
experiments 0, 1^], which study resonantly interacting 
ultracold ^Li atoms in a trap with a population imbal- 
ance of different spin components. The pairing super- 
fluidity in the Fermi gases with population imbalance 
between the two spin components have been studied be- 
fore in different physical contexts, m ostly in the weakly 
interacting regime 0, IE S IE IE HQ- Some exotic 
phases have been proposed to exist due to the compe- 
tition between the Cooper pairing and the population 
imbalance |a, S S U^ ■ The recent experiments near a 
Feshbach resonance have raised strong theoretical inter- 
est in studying the phase configuration of such a system 
under a potential trap in the strongly interacting region 

[iiiiiiaiiiiiiiaiiiiii. 

In this work, we map out the detailed phase diagram 
for fcrmionic atoms in a trap with population imbalance, 
both at zero and at finite temperature. We use the same 
theoretical method as proposed in Ref. 13], which cor- 
responds to a generalization of the self-consistent GqG 
diagram scheme 19] from the equal population case to 
the case with population imbalance. At zero tempera- 
ture, this method reduces to the mean-field approxima- 
tion for the crossover theory 13, 15, 16, 20]; while at fi- 
nite temperature, it includes a pseudogap in addition to 
the superfluid order parameter. To avoid subtle unstable 
solutions for the ground state of this system, we directly 
minimize the thermodynamic potential under the local 
density approximation instead of using the gap equation. 

One of the difficulties to map out the detailed phase 
diagram lies in the fact that the properties of the system 
seem to depend on so many different parameters. For 
instance, we expect in general, several different phases to 
exist from the trap center to the edge, with their char- 
acters and boundaries determined by the temperature of 
the system, the population imbalance, the magnetic field 
detuning, the atom specie, the trap frequencies along the 
three spatial directions, and the total atom number. It is 
difficult to calculate the phase distribution for all these 
different parameters. Fortunately, similar to the homoge- 



neous system with equal spin populations, there exists a 
nice universality for this more involved system. Although 
the properties of the system depend on all the parameters 
mentioned above, the dependence is through some dimen- 
sionless combinations of the physical parameters. As a 
result, the phase diagram is completely fixed by three 
dimensionless parameters after rescaling: the normalized 
temperature, the dimensionless interaction strength, and 
the population imbalance. In particular, at zero tem- 
perature and at the resonance point, the phase diagram 
only depends on a single parameter: the population im- 
balance. In this universal frame, the variations in the 
trap (anisotropic in general) or in the atom number do 
not induce any further complexity for the description of 
the system. 

To fix the phase diagram, we calculate under various 
interaction strengths and temperatures, the distribution 
of the system's phases from the trap center to the edge as 
a function of the population imbalance. The main results 
are shown in Fig. 1 and 3. In the following, we first 
give the universal equations of state in Sec. II, written 
in terms of the three dimensionless parameters. In Sec. 
Ill and IV, we present our main calculation results with 
detailed discussions. 



II. THE FORMALISM FOR A TRAPPED 

FERMI GAS WITH POPULATION IMBALANCE 

AND THE UNIVERSALITY 

As the population of the closed channel molecules is ex- 
ceedingly small near a wide Feshbach resonance |l3, ^M ' 
it is sufficient to use the following single-channel Hamil- 
tonian to describe the Fermi gas in the near resonance 
region: 



H = 



k.cr 



(ek - /"<j)ak.<jak,, 



(1) 



(U/V) 



E 

q,k,k' 



"q/2+k,T q/2- 



.k^xaq/2-k'aaq/2+k',T 



where ek = k"^ /{2m) (to is the atom mass and h — 1), 
/i(j is the chemical potential for the spin-cr component 
(cr =T, i labels the two spin states), V is the quantization 
volume, Oj^ ^ is the creation operator for the fermionic 
atoms. The bare atom-atom interaction rate U is con- 
nected with the physical one Up — Anas/m (a^ is the 



atomic scattering length) through the standard renormal- 
ization relation l/U = l/[/p- (1/V) ^^ 1/ (2ek) ilSi- We 
take the local density approximation so that /X| = jir + h, 
Mi = Mr — ^j Mr = M ~ ^ (r), where V (r) is the exter- 
nal trap potential (slowly varying in r). The chemical 
potential /i at the trap center and the chemical poten- 
tial imbalance h are determined from the total atom 
number N = N-^ + N^ and the population imbalance 
/3 = |iV| — Ni\ /N through the ntimber equations below. 

As has been shown in Ref. [13, under the local den- 
sity approximation, the thermodynamic potential fJ = 
— Tln[tr(e^^/-'")] corresponding to Hamiltonian (1) has 
the following expression: 



^IV = -|A|Vt/p-(r/V)^{ln[l + exp(-|i?kil/r)] 

k 

+ ln[l + cxp(- li^kT I m - I A| V(2ekT) (2) 
-[ek - Mi - 0{E^i)E^i + 0(-Ski)i?ki]/T}, 



where the gap A at zero temperature is given by the 
order parameter As = C/X]k('^-k,iak,T); and at finite 
temperature should be understood as the total gap, with 

|A| = 



is I -|- I Apg I , where Apg is the pseudogap com- 
, from the contribution of non-condensed Cooper pairs 
The 6'-function is defined as 9 [x) ~ 1 for a; > 
and 9 (x) = otherwise. Without loss of generality, we 
take /i > so that N-\^ > N^ always. Note that different 
from the equal-population case, the quasi-particle excita- 
tion energies E^ko- are different for the a — '\, I branches: 



^kT,i = Y (fk — Mr)^ + |A| =F h. In the case of /i > 0, 
E'ki is always positive; while there exists certain param- 
eter regions where the sign of £^kT becomes momentum 
dependent and is negative in the range fc_ < |k| < fc+, 

where fc^ = max 0, 2m(/ir ± \/h^--\A^) . In this mo- 
mentum range, the atoms stay unpaired, which corre- 
sponds to the so-called breached pair state |3, 123 • In 
deriving the thermodynamic potential (2), we have ne- 
glected the non-zero-momentum pairing (the FFLO state 
[5i|, with the pair momentum q ^^0). This is motivated 
by the fact that the FFLO state is stable only within a 
narrow parameter window in the BCS regio njlft Il3. l2^| 
and is absent in the recent ^Li experiments 0, 01 ■ 

From the thermodynamic potential, one can get the 
gap equation from the condition dfl/dA = 0. However, 
as discussed in Ref. 1^ , in the presence of a population 
imbalance, the thermodynamic potential has a double 
well structure, and the gap equation may give unstable 
solutions. Therefore, instead of solving the gap equation, 
we directly minimize the thermodynamic potential Q to 
find its global minimum with respect to the gap param- 
eter A. To fix the chemical potentials fia in Eq. (2), we 
need to use the number equations, derived from the rela- 
tions dn/diia- — — "■ro-V, where Ura- is the number density 
of the spin-cr component at the position r. The number 



equations have the form 

nr. = ^ E["k/(^k..) + «^/(-i?k,-.)], (3) 

k 

where the parameters u^. — (_Ek + (fk — Mr))/2-Ek, v^ = 

(Ek - (ek - /ir))/2Ek, Ek = Y^(^i~^O^T^, the 
Fermi distribution f{E) = 1/ (l -f e^/"^), and for con- 
venience, we take — 1=1 and vice versa. The atom den- 
sities nr| and n^i arc connected with the total atom 
number and the population imbalance through N — 
J d^r{nr^ + Uri), and/? = J d^rSur/N {Srij. = nr|— Urx). 
Given an external potential V (r) , we can determine 
the properties of the system from Eqs. (2) and (3). 
However, the solution in general depends on too many 
external parameters. If the trapping potential is har- 
monic with the form ^(r) = X^i 5'^'^?^? (^ ~ 2;,j/,z, 
anisotropic in general) , the properties of the system ob- 
viously will depend on the temperature T, the effec- 
tive scattering length Os, the atom mass to, the total 
atom number TV, the population imbalance /3, and the 
trap frequencies uji along the three spatial dimensions. 
This much involved dependence can be significantly sim- 
plified if we transform the set of equations above into 
the dimensionless form. For that purpose, we choose 
the unit of energy to be the Fermi energy (Ep) at the 
center of the trap for N non-interacting fermions with 
equal population for the two spin components. Under 
the local density approximation, one can easily figure out 
Ep = (SNuJxUjyLUz)'^ from its definition. The harmonic 
trap in the unit of Ep can be expressed in the standard 
dimensionless form V{r)/Ep = X^i'^j where r^ = ri/Ri, 
with the Thomas-Fermi radius Ri = \/2Ep/mujf along 
the ith direction. With these, the number equations in 
(3) are cast into the following dimensionless form: 

l±P=^ j d^Td^k[ulf{E^^„) + vlf{-E^,_^)], (4) 

where the energies and the momenta are in the units of 
Ep and kp [kp = J2mEp/fi ), respectively. The di- 
mensionless thermodynamic potential has the same form 
as Eq. (2), except that the interaction strength Up is 
replaced by the dimensionless one -^kpUs, and that all 
the energies are normalized by the chosen unit Ep (cor- 
respondingly, k by kp and T by Tp = Ep/kB)- 

From these dimensionless equations, it becomes obvi- 
ous that the properties of the system depend only on the 
three dimensionless parameters T jTp^ kpUg and /3. The 
system will have the same properties as long as these 
three parameters are the same, whether it is for differ- 
ent atom species, or with different total atom numbers, 
or in traps with different trapping frequencies (uJi). This 
shows that the properties of this inhomogeneous system 
still have nice universality, similar to the case of a ho- 
mogeneous Fermi gas without the population imbalance 
|20l |23| , where the system can be characterized by two 



dimensionless parameters kpa^ and T/Tp. Comparing 
with the case of a homogeneous gas, we see that in the 
current case, the presence of (anisotropic) traps with var- 
ious atom numbers do not add complexity to the descrip- 
tion of the system. This nice feature comes from the local 
density approximation and the assumption of harmonic 
traps 245 , and is independent of the particular approxi- 
mation schemes in deriving the equations of state. 



III. PHASE BOUNDARIES AT ZERO 
TEMPERATURE 



Following the formalism outlined in the previous sec- 
tion, we first map out the phase boundaries for trapped 
fcrmions at zero temperature. The different phases in 
the trap can be identified from the gap A and the chem- 
ical potentials /i^ of the two spin components. From the 
solutions of A and /ig- , one can immediately get the den- 
sity profiles of the trapped atoms. There are in general 
four possible phases in this system [l^ [I^ : (i) a super- 
fluid (SF) state with A > and no Fermi surface in the 
momentum space; (ii) a breached pair (BPl) state with 
A > 0, one Fermi surface in the momentum space, and 
gapless fermionic excitations J25| ; (iii) a normal polarized 
(NP) state with A = and one Fermi surface (/i| > 0, 
Hi < 0); and (iv) a normal mixed (NM) state with A = 
and two Fermi surfaces (/X| > /i| > 0). At zero tempera- 
ture, the density profiles of the two spin components are 
identical in the SF phase, but are different in the BPl 
phase. This can serve as the measure to distinguish these 
two phases at T = 0. At finite T, both of the phases have 
finite polarization (excess fermions) with different density 
profiles for the two spin components [l3|- In this case, 
they can only be distinguished by the existence of the 
Fermi surface in the BPl phase (see the note (251). 

Figure 1 shows the detailed phase distribution in the 
trap at different interaction strengths {kpas)~^ (corre- 
sponding to different magnetic field detunings). On the 
BEC side of the resonance (Fig. 1(a), with {kpag)^^ = 3) 
and at small but nonzero population imbalances, the 
Fermi gas separates into three phases in the trap: an 
SF core at the center, a BPl phase in the middle, and 
an NP phase at the edge of the trap with only the ma- 
jority spin component. As the imbalance parameter (3 
increases, the superfluid core becomes smaller until it 
vanishes at a critical imbalance, beyond which only the 
BPl phase and the normal phase exist in the trap. This 
critical imbalance where the phase transition from the 
SF phase to the BPl phase occurs becomes greater to- 
wards the resonance, while the parameter range of the 
BPl phase shrinks (Fig. l(b)(c)). ki {kpa^y^ = I, the 
BPl phase only exists at small population imbalances in 
a slim region bordering the SF and the NP phase (Fig. 
1(c)). Within our numerical resolution {5(3 ^ ±1%), the 
BPl phase disappears at roughly [kpas)^^ ^ 0.5, on the 
BEC side of the resonance. Note that from numerical 
analysis, we find that the unpaired fermions in this BPl 



phase are within the momentum range < |k| < A:_|_ 
(which implies |/ir| < ^h"^ — |Ap in this region), where 
fc-i_ is given in the previous section. In this BPl phase, 
we therefore have the following phase separation picture 
in the momentum space: the paired fermions in the su- 
perfluid fill the outside momentum shell, while the un- 
paired fermions of the majority spin component occupy 
the states inside the Fermi ball with |k| < /c+. 

Once the BPl phase disappears, the trap is left with 
only the SF phase and the NP phase over a certain range 
of the interaction strength (Fig. 1(d), with {kpas)~^ ~ 
0.5). This continues till a new normal state, the NM 
phase shows up at large population imbalances at roughly 
{kpas)~^ ^ 0.1 (Fig. 1(e)), where fermions of both spin 
components show up in the normal gas, with different 
Fermi surfaces for different spins. The range of the NM 
phase grows towards the resonance, and at resonance or 
on the BCS side (Fig. l(f)(g)(h)), the Fermi gas typi- 
cally phase separates into three regions at small /3: the 
SF phase, the NM phase and the NP phase. The super- 
fluid phase disappears at a critical population imbalance, 
where the gas undergoes a phase transition from the su- 
perfluid phase to the normal state. Qualitatively, this 
picture agrees pretty well with the recent experimental 
findings ,2j , although quantitatively, the mean-field type 
approximation at T = may somewhat overestimate the 
critical population imbalance for the disappearance of the 
SF core at the trap center. 

To give more detailed information of this system, we 
also show in Fig. 2 the chemical potential jjlq at the 
trap center and the chemical potential difference h as 
functions of the population imbalance /3 at the corre- 
sponding interaction strengths {kpas)~^ . Note that the 
chemical potential difference h between the two spin com- 
ponents does not change across the trap. The local chem- 
ical potential /ir changes, but given /io and the trap 
potential V{y), it changes through the simple relation 
A'r = Mo — ^(r) under the local density approximation. 
With the information given in Fig. 2, we know the local 
chemical potential /iro- for each spin component cr, and it 
becomes straightforward to calculate other properties of 
the system, such as the density profiles. 



IV. PHASE BOUNDARIES AT FINITE 
TEMPERATURE 

At finite temperature, we calculate the total gap (A) 
in the trap, from which we get the boundary between the 
gapped region and the gapless region. As the order pa- 
rameter for the superfluid phase A^ is smaller than the 
total gap at finite temperature, the boundary between 
the gapped and the gapless regions serves as an upper 
bound for the superfluid phase |2y|. At flnite tempera- 
ture, due to the thermal excitations, the atom density 
profiles change more smoothly as one goes from the trap 
center to the edge, so one cannot easily use the disconti- 
nuity of the density profiles to fix the boundary between 
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FIG. 1: (Color Online) The zero temperature phase boundaries of a polarized Fermi gas in a trap with (a)(A;Fas) ^ = 3, (b) 
{kpas)-^ = 1.5, (c) {kpas)-^ = 1, (d) {kpas)-^ = 0.5, (e) {kpas)-^ = 0.09, (f) {kpas)-^ = 0, (g) [kpas)-^ = -0.4, (h) 
{kpcLs)'^ ~ —1. The solid lines mark the phase boundaries between the superfluid region (SF/BPl) and the normal region 
(NM/NP); the dashed lines are the phase boundaries between the SF phase and the BPl phase; the dotted lines show the range 
of the minority spin component in the normal phase, which are effectively the phase boundaries between the NM and the NP 
phase; and the dash-dotted lines mark the range of the majority spin component in the normal phase. The trap radius r is in 
the units of the Thomas-Fermi radius for the corresponding directions. 















(a) 






^ 








\ 
















~~- 


-h 


- 


- 


- 


(kpar'=3 



1.8 




















. 


-2 










^ 












2.2 




















- 


2.4 


\ 


















- 


2.6 
28 
-3 

3.2 


\ 


^^ 


- 


- 


-h 


- 


- 


- 


— 





3.4 




















{kpar'=1.5 . 



0.3 1 



0.4 0.6 O.S 




0.8 
0.7 
0.6 
OJ 
0.4 
03 
0.2 
0.1 




















(e) 






\ 
























































^ 


" 


— 










^ 


^ — 














^ 


h 










" 


/ 














. 


/ 




























<v, 


"'=0.09 















(<i) 








^ - 


-- 


- " 


-' 




-- 


h 










/ 














/ 














? 












(k^a)- 


=0.5 " 






\ 



■ 
















(f) 


0.8 


L- 





















0.7 














\ 


^ 


0.6 
















- 


O.i 














^ 


- " 


0.4 














^ — " 


" 


0.3 








^ 


- 


- 




- 


0.2 




/ 


/ 








h 


" 


0.1 


'1 














(k a )"'=0 



0.8 1 



0.8 1 























(g) 


0.8 






\ 














0.6 




















/ 


0.4 

02 


/ 


___ ^ - 


h 


- 


- 


- 


- 


- 


k^a 


r'=-o.4 













"^ ^~^^^^_^ ,h) 


0.8 










\ 


0.6 










J 


0.4 

0.2 


^ 







— 





0.8 1 



0.2 0.4 



0.8 1 



FIG. 2: (Color Online) The chemical potential /xo and the chemical potential difference h as functions of population imbalance 
at zero temperature for: (a)(A:Fas)~^ = 3, (b) (feFfls)"^ = 1-5, (c) (feFfls)"^ = 1, (d) (/cFfls)"^ = 0.5, (e) (feFffls)"^ ~ 0.09, (f) 
(^Ffls)"^ = 0, (g) (/sFfls)^^ = —0.4, (h) (fcFfls)"^ ~ —1- The solid lines represent ^,q/Ef^ the average chemical potential at 
the center of the trap; and the dashed lines are h/ Ef { — h/Ep in (a-c) for better comparison with jio/Ep), where the unit of 
energy Ef is defined in Sec. II. 
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FIG. 3: (Color Online) The finite temperature phase boundaries of a polarized fermion gas in a trap with (a)(fcFas) ^ = 3, 
(b) (fcFfl,)-' = 1.5, (c) (kpas)-^ = 1, (d) (kpas)-'' = 0.5, (e) (ftFa,)"' = 0.1, (f) (/cFa,)"' = 0, (g) (fcFa,)"' = -0.4, (h) 
{k — Fas)"^ = —1. The solid lines mark the phase boundaries between the gapped and the gapless region; the dashed lines 
are the phase boundaries between the gapped phases without a Fermi surface and the ones that have; the dotted lines show 
the zero point of chemical potential of the minority spin component in the normal phase; and the dash-dotted lines show the 
zero point of the chemical potential of the majority spin component in the normal phase. The temperature is taken to be 
T = 0.22Tf, which corresponds to a real temperature T ~ SOOnK for 2.7 x 10^ fermionic ^Li atoms in a cigar shaped trap 
with u!z ~ 23Hz,uJx = ojy ^ llOHz Q- 
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FIG. 4: (Color Online) Same as Fig. 2, but [i, and h are shown as functions of population imbalance at finite temperature 
(T = 0.22Tf) with: (a)(fcFa,)~^ = 3, (b) {kpas)'^ = 1.5, (c) [kpasY^ = 1, (d) (kpas)'^ = 0.5, (e) {kpasY^ = 0.1, (f) 
(fcFds)"^ = 0, (g) (fcFffls)"^ = —0.4, (h) (fc — FasY^ — —1. The solid lines represent ^,/ Ef\ and the dashed lines are h/Ep 
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the BPl phase and the SF phase, nor the boundary be- 
tween the different normal states (NM/NP). However, 
one can stih look at the changes of the Fermi surfaces 
in the momentum space, which provides an unambigu- 
ous signal to determine the boundaries between different 
phase/regions. We should caution, however, that these 
boundaries do not necessarily correspond to sharp edges 
in the atom density profiles. For instance, in the normal 
phase with A = 0, when one of the chemical potentials, 
say /ir|, changes its sign at the trap edge, which im- 
plies the disappearance of the Fermi surface for the spin 
up atoms as one moves outside (under the local density 
approximation), the density distribution of the spin up 
atoms does not vanish at that point. Instead, it will fol- 
low an exponential decay as the corresponding chemical 
potential turns negative. So, although these boundaries 
do not represent sharp spatial edges of the corresponding 
spin components, they indeed give a good estimation of 
the ranges. 

The results of our calculation at finite temperature are 
shown in Fig. 3. We take the temperature T = 0.22Tp, 
which corresponds to T ~ 300nK with the experimen- 
tal parameters in Ref. U (see caption of Fig. 3). The 
phase diagrams are qualitatively similar to those in Fig. 
1, but there are several important differences at finite 
T which deserve to be emphasized. First of all, on 
the BEC side, the finite temperature BPl region does 
not show up until a critical population imbalance (see 
Fig. 3(a-d)). The basic reason behind this feature is 
that at finite temperature, the population imbalance can 
be carried by the quasiparticle excitations in the con- 
ventional BCS superfluid state, as we have emphasized 
in |l3l |. Therefore, different from the zero temperature 
case, the conventional BCS state becomes partially po- 
larized when there is a population imbalance, which sig- 
nificantly relaxes the competition between the population 
imbalance and the Cooper pairing, and which makes the 
momentum-space phase separation between the paired 
state and the excess fermions (the BPl phase) unneces- 
sary at small population imbalances. Secondly, the range 
of the BPl region changes significantly from the zero tem- 
perature case. At T = 0.22Tp, the BPl region disappears 
near (fc^^as)"^ ~ 0.15, where there already exists an NM 
region in the phase diagram (see Fig. 3(d) (e)), which 
initially appears near (fci^Og)"^ ~ 0.3 for small popula- 
tion imbalance. Last but not least, at resonance or on 



the BCS side, the critical population imbalance at which 
the gapped region disappears at the center of the trap, 
which is an upper bound for the actual critical imbalance 
bordering the SF phase and the non-condensed pairs, be- 
comes significantly smaller than the critical population 
imbalance at zero temperature (see Fig. 3(f)(g)(h)). The 
results above demonstrate that the phase boundaries, as 
well as the evolution of the different phases as the field 
sweeps across the crossover region, can be significantly 
affected by temperature. As all the experiments are nec- 
essarily done at a finite temperature, this suggests that 
it may be important to take the thermal effects into ac- 
count for a quantitative interpretation of the experimen- 
tal data. 

At finite temperature, we also calculate the chemical 
potentials /iro- as functions of the population imbalance 
/3 at various interaction strengths (fci^Os)"^, in order to 
provide detailed information of this system. The results 
are shown in Fig. 4. The properties can be directly 
read from the figure, and as the main features there are 
qualitatively similar to those in Fig. 2 for the zero tem- 
perature case, we neglect the detailed discussion here. 



V. SUMMARY 

In summary, we provide detailed calculations of the 
phase diagrams for a trapped Fermi gas with population 
imbalance over the entire BCS-BEC crossover region at 
zero and at finite temperature. The calculation is done 
with the self-consistent GqG diagram scheme. We also 
emphasize the importance of using the universal dimen- 
sionless equations. The properties of the system then 
only depend on three dimensionless parameters, from 
which we can calculate the universal phase diagrams valid 
for different atom species with different atom numbers 
under various trap configurations. The main results of 
our calculation are shown in Fig. 1 and 3, with their 
prominent features discussed in detail in the correspond- 
ing sections. 
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If the trap has anharmonicity which is a constant ex- 
pressed in the unit of the Thomas- Fermi radius (see M. 
W. Zwierlein, W. Ketterle, lcond-m at/0603489 ), we will 
still have the same universality as discussed in the text. 
Note that this state is also called the polarized (mag- 
netized) superfluid state in the literature 10]. However, 
at finite temperature, the normal BCS superfluid state 
also has finite polarization as emphasized in Ref. p^ . 
To distinguish between these different kinds of super- 
fluid phases, we call the non-BCS superfluid as the BPl 
phase (abbreviation for the breached pair phase with 
one Fermi surface), similar to the terminology of Ref. 
[J. The breached pair state was studied in Ref. Q on 
the BCS side, where it has two coexisting Fermi surfaces 
(BP2 phase). Other than that difference, the properties 
of the BPl and BP2 phases are pretty much the same. 
For instance, they are both superfluid states with gap- 
less fermionic excitations, and are both homogeneous in 
the real space while phase separated in the momentum 
space. 

To exactly flx the phase boundary for the superfluid 
phase at flnite T, one needs to derive the dispersion re- 
lation for the pair excitations to split A into As and 
Apg. This dispersion relation is known for the equal- 
population case 19] , but has not yet been derived for the 
case with population imbalance within the self-consistent 
GqG scheme. However, for the calculation of the total gap 
A and the atom density proflles, one does not need to flx 
that dispersion relation [l^. 



